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Abstract 

It is now widely understood that all animals engage in complex interactions with bacteria (or microbes) throughout their 
various life stages. This ancient exchange can involve cooperation and has resulted in a wide range of evolved host- 
microbial interdependencies, including those observed in the gut. Ciona intestinalis, a filter-feeding basal chordate and 
classic developmental model that can be experimentally manipulated, is being employed to help define these relationships. 
Ciona larvae are first exposed internally to microbes upon the initiation of feeding in metamorphosed individuals; however, 
whether or not these microbes subsequently colonize the gut and whether or not Ciona forms relationships with specific 
bacteria in the gut remains unknown. In this report, we show that the Ciona gut not only is colonized by a complex 
community of bacteria, but also that samples from three geographically isolated populations reveal striking similarity in 
abundant operational taxonomic units (OTUs) consistent with the selection of a core community by the gut ecosystem. 
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Introduction 

Sequential colonization of newborns has ancient origins and has 
resulted in the evolution of host-symbiont interdependencies [1- 
3] . Studies comparing wildtype animals to those maintained under 
sterile conditions reveal bacterial influences on host development, 
tissue maturation and metabolic capabilities, indicating complex 
functional associations between host and microbe [4] that may 
depend on a continuous exposure to microbes [5-7]. Recent 
studies indicate that not all animals are dependent on a core 
microbiota. Instead, some maintain residents that fulfill core 
requirements (i.e., a functional core), suggesting that some 
interactions or interdependencies with bacteria may be less 
species-specific than considered previously [8,9]. In either case, 
it appears that most animals can be considered metaorganisms 
that rely on a multitude of microbial associations [3,4,10-12]; even 
in marine sponges, which lack distinct body compartments, 
species-specific microbiota appear to be maintained [3,13]. Ciona 
intestinalis, a sessile filter-feeding marine invertebrate chordate, 
possess mechanisms that likely help it regulate host-microbe 
interactions in the gut [14,15]. In order to better understand and 



describe the composition of bacterial communities associated with 
the Ciona gut, we sampled animals from three geographically 
disparate populations and undertook several different approaches 
to define their respective gut microbiota. 

Materials and Methods 

Ethics statement 

The research described herein was performed on the marine 
invertebrate, Ciona intestinalis and did not involve human or other 
primate subjects or mammals of any form. Ciona is not protected 
by any environmental agency in the United States (US) or in Italy. 
The collection services contracted in this study maintain current 
permits and licenses for collection and distribution of marine 
invertebrates to academic institutions; special permission was not 
required to collect these animals. Handling of live animals was in 
accordance with the guidelines of our academic institutions. 
Animals were recovered and brought to the laboratory alive and 
maintained in clean water with aeration. DNA was collected from 
the gut of dissected animals. Animal waste products were disposed 
of appropriately. 
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Animal Collection 

Ciona intestinalis were acquired in April 2012 from three distinct 
populations: near Mission Bay, San Diego (S. LePage, M-Rep, 
Carlsbad, CA), at the Cape Cod Canal Basin near Woods Hole 
(Animal Collection Services at Marine Biological Laboratories, 
MA), and from Fusaro Lake in Naples, Italy (Thyrrenian Sea, 
Italy, Animal Collection Services, Stazione Zoologica Anton 
Dohrn). The starved animals are designated WH(a), SD(a), and 
N(a) and two sets of unstarved individual animals, which were 
collected previously in April 201 1 from the same location as WH 
and SD and are designated here as WH(b),(c) and SD(b),(c) (gut 
contents reflect diet at collection site), also are included in the 
analysis. 

All animals were sampled immediately (at random) from the 
population, bagged and shipped overnight to our laboratory (in 
the US or Italy), and processed immediately for 16S rRNA gene 
PCR. Animals to be starved were maintained in 0.2 micron 
filtered seawater (FSW) at arrival, see below. 

Animal and sample handling 

Starved animals (4—5 animals, pooled for DNA) were main- 
tained in sterile containers at room temperature in 0.2 micron 
FSW (with complete water changes every 4—6 hours) for 72- 
96 hours to void gut content. The gut (stomach and intestine) from 
each animal was dissected, homogenized and total DNA was 
isolated (Qiagen Stool Kit) using aseptic technique and benchtop 
clean units. The starved samples ("a" samples), one pooled set 
from each of the three populations (labeled internally as Gi- 121, 
Cil22, and Cil26 or WH, SD, and Naples, respectively), were 
sequenced separately as described below and used as templates to 
characterize the bacterial community structure. Unstarved 
animals, Ci-41, -42 (WH) and Ci-47, 48 (SD) were processed in 
the same manner, and designated in the text as samples "b" and 
"c," respectively, for each population. 

16SrRNA PCR Clone Library and Sequencing 

PCR products were generated with 27F and 1492 primers [16- 
1 9] using Promega 2 x Master Mix (Promega, Madison, WI, USA) 
under aseptic conditions. Products were cloned into TOPO-based 
cloning vectors (Invitrogen/Life Technologies, Carlsbad, CA, 
USA). Bacterial clones were screened for inserts by colony-PCR 
using vector-specific primers and sequenced on an ABI3730 
automated sequencer (Applied Biosystems/Life Technologies, 
Carlsbad, CA, USA). 

16SrRNA Library Preparation and Amplicon Sequencing 
Method 

DNA template was amplified using the 515F/806R region of 
the 16S rRNA gene using primers and cycling conditions modified 
slightly from those outlined in Caporaso et al. 201 1 [20], adapted 
for the Illumina HiSeq2000 and MiSeq by adding nine extra bases 
in the adapter region of the forward amplification primer that 
support paired-end sequencing. Briefly, the V4-5 region of the 16S 
rRNA gene was amplified with region-specific primers that 
included the Illumina flowcell adapter sequences. The reverse 
amplification primer also contained a twelve base barcode 
sequence that supports pooling of up to 2,167 different samples 
in each lane. PCR reactions were conducted in triplicate and 
products were pooled. Each pool was then quantified using 
PicoGreen (Invitrogen) and a plate reader. Once quantified, 
different volumes of each of the products were pooled into a single 
tube so the pool consists of an equal amount of each amplicon and 
cleaned using the UltraClean® PCR Clean-Up Kit (MoBIO, 



Carlsbad, CA, USA). Amplicons then were sequenced in a 
151 bpxl2 bp x 151 bp MiSeq run employing custom sequencing 
primers and procedures described in the supplementary methods 
of Caporaso et al. 2012 [21]. 

Sequence Analysis 

All sequence analyses were completed in QJIME 1.6.0-dev 
[22,23]. OTUs were selected based on 97% sequence similarity 
and taxonomic data were assigned to each representative sequence 
through the classification algorithm of the Ribosomal Database 
Project (RDP). Chloroplast sequences derived from eukaryotic 
endosymbionts were discarded, and sample libraries were rarified 
to a common depth of 40,528 reads before comparison. All cluster 
analyses are based on weighted unifrac distance because all 
samples were found to be dominated by a small number of high 
abundance OTUs; the extent to which those OTUs dominate the 
samples therefore is more relevant than weighting low abundance 
taxa disproportionality by comparing samples based on the 
presence or absence of OTUs. 

Illumina-derived core OTUs and clone library, Sanger- 
sequenced, OTUs were aligned using Geneious version 6.1.6 
(available from Biomatters, www.geneious.com). OTUs were 
compared by phylogenetic reconstruction using the Minimum 
Evolution method in Mega 5.2.2 [24] with 1000 bootstrap 
replicates and the Maximum Composite Likelihood nucleotide 
substitution model, to include both transitions and transversions of 
non-protein coding DNA sequences. Uniform rates were main- 
tained, gaps were treated by complete deletion, and trees were 
inferred by Close-Neighbor-Interchange. 

Results and Discussion 

Evidence for core OTUs in the Ciona gut 

Each sample was sequenced to a depth of 40,528 reads (data is 
deposited in MG-RAST under the following IDs: 4539627.3, 
4539628.3, 4539629.3, 4539630.3, 4539624.3, 4539625.3, 
4539626.3), which collectively encompassed over 9,700 unique 
OTUs dominated by 8 bacterial phyla (Figure 1A, ID). Despite 
this high degree of community complexity, the 1 2 most abundant 
OTUs in the study were sufficient to capture at least half of the 
microbial diversity in all but one sample, suggesting the presence 
of a core gut community (Figure IB). These 12 OTUs, 9 of which 
were detected in all 7 samples and all of which were detected in at 
least 5, each contained at least 2% of all reads in the study. 
Although a relatively small percentage of all detected OTUs were 
found in both starved and unstarved samples, those OTUs 
comprised the vast majority of observed 16S reads, especially after 
discounting singleton OTUs (OTUs containing only a single 
sequence across all samples) (Figure 1C). Accordingly, observed 
differences in community composition are largely attributable to 
very low abundance OTUs, and may be an artifact of limited 
sampling depth. Rarefaction curves depict slightly greater 
observed diversity in the samples taken from unstarved Ciona 
relative to their starved counterparts, although the difference is not 
statistically significant (Figure ID). Overall, we found 35 OTUs 
shared among all 7 samples and 200 OTUs present in a least 75% 
(5) of the samples (Table 1 and Table SI). 

The gut microbiota of Ciona is predominated by Gram-negative 
bacteria and among the core OTUs (Table 1 ), Proteobacteria were 
the most abundant (28 OTUs or 80%), followed by Bacteroidetes 
(4 OTUs), Fusobacteria (2 OTUs), and Verrucomicrobia (1 
OTU). Together, this assemblage represents 13 families of bacteria 
(Table 2), likely conferring broad metabolic potential to the host 
gut [25,26]. In addition, OTUs related to the Gram-positive 
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Figure 1. Taxonomic summary of Ciona samples by lllumina sequencing of 16S rRNA (alpha diversity). (A) Relative abundances of 
sequences classified to phylum, with Proteobacteria split by class. Only phyla containing at least 1% of reads across all samples are designated by 
color. (B) Relative abundances of the 12 OTUs containing at least 2% of all reads across samples. OTUs are classified to the finest taxonomic level 
possible, with class indicated in parentheses where appropriate. The dendrograms at left of the figures depict UPGMA clustering of the seven 
samples based on weighted community similarity. (C) Distribution of OTUs and their assigned sequences between starved and unstarved samples. 
Although a relatively small percent of all detected OTUs were found in both starved and unstarved samples, those OTUs comprised the vast majority 
of observed 16S reads, especially after discounting singleton OTUs. (D) Rarefaction curves for starved and unstarved communities, depicting the 
greater number of observed species in samples which were not starved than in those that were. Error bars are standard deviation. 
doi:1 0.1 371 /journal.pone.0093386.g001 



Firmicutes and Actinobacteria are likely important members of the 
core community and contribute several abundant OTUs in more 
than 75% of the samples. As evidenced by UPGMA clustering of 
community similarity (Figure 1 A) and network clustering of sample 
OTUs (data not shown), there is no observed correlation either 
between the geographic sampling location or the sampling method 
(starvation) and the beta diversity relationships between samples 
(Table S2) as indicated by the G-test of independence producing a 



p value of 1 [23]. This G-test result indicates that despite whatever 
apparent clustering exists, two samples from the same location are 
not more likely to share an OTU than are samples from different 
locations. Similar findings were noted for starved versus unstarved 
animals. The detection of identical OTUs in the guts of Ciona from 
three geographically disparate populations is consistent with a 
global marine microbial seedbank [2 7] .Taken together, these 
observations support the notion of a core microbial community in 
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Table 1. Core OTUs (35) found in all samples* from three geographically disparate populations. 





OTU ID 


N(a) 


WH(a) 


SD(a) 


WH(b) 


SD(b) 


WH(c) 


SD(c) 


% Total 


Phylum 


Family 


8443 


5219 


965 


9503 


1467 


5028 


3545 


21687 


16.7% 


Proteobacteria 


Endozoicimonaceae 


8554 


2 


13633 


4939 


12420 


81 


124 


38 


1 1 .0% 


Fusobacteria 


Fusobacteriaceae 


8942 


14 


9300 


1417 


1685 


38 


1388 


23 


4.9% 


Proteobacteria 


Vibrionaceae 


6994 


8147 


7 


404 


451 


872 


199 


2402 


4.4% 


Proteobacteria 


Unresolved (Alteromonadales) 


10014 


2 


28 


196 


5974 


3 


5630 


1 


4.2% 


Proteobacteria 


Psychromonadaceae 


1955 


3043 


1022 


2161 


906 


680 


38 


1840 


3.4% 


Fusobacteria 


Fusobacteriaceae 


717 


8324 


2 


27 


64 


23 


59 


133 


3.0% 


Proteobacteria 


Unresolved (Gamma-) 


2967 


1396 


468 


1496 


1219 


324 


735 


1271 


2.4% 


Proteobacteria 


Vibrionaceae 


4938 


34 


174 


1450 


1684 


352 


156 


874 


1.7% 


Proteobacteria 


Desulfobulbaceae 


3583 


5 


3 


40 


1 


3180 


12 


50 


1 .2% 


Proteobacteria 


Unresolved (Alpha-) 


6378 


1 


26 


1286 


172 


250 


1413 


84 


1.1% 


Proteobacteria 


Campylobacteraceae 


9665 


2 


1864 


276 


596 


16 


269 


20 


1.1% 


Proteobacteria 


Shewanellaceae 


9259 


5 


898 


81 


245 


665 


172 


92 


0.8% 


Bacteroidetes 


Unresolved 


234 


1 


572 


290 


2 


873 


115 


81 


0.7% 


Bacteroidetes 


Flavobacteriaceae 


5491 


1 


257 


1242 


2 


53 


141 


37 


0.6% 


Proteobacteria 


Oceanospirillaceae 


1617 


293 


327 


79 


239 


114 


257 


49 


0.5% 


Proteobacteria 


Vibrionaceae 


838 


253 


17 


86 


74 


96 


277 


30 


0.3% 


Bacteroidetes 


Flavobacteriaceae 


1263 


20 


27 


8 


2 


592 


16 


40 


0.2% 


Proteobacteria 


Rhodobacteraceae 


2485 


12 


14 


67 


6 


467 


58 


68 


0.2% 


Proteobacteria 


Rhodobacteraceae 


4282 


17 


265 


34 


4 


59 


230 


63 


0.2% 


Proteobacteria 


Pseudoalteromonadaceae 


1329 


1 


5 


83 


2 


241 


1 


95 


0.2% 


Bacteroidetes 


Unresolved 


8677 


3 


2 


3 


12 


295 


57 


33 


0.1% 


Proteobacteria 


Rhodobacteraceae 


6190 


2 


6 


95 


83 


46 


163 


9 


0.1% 


Bacteroidetes 


Campylobacteraceae 


1489 


5 


12 


98 


15 


34 


32 


39 


0.1% 


Proteobacteria 


Alteromonadaceae 


6952 


11 


7 


38 


4 


108 


20 


31 


0.1% 


Proteobacteria 


Rhodobacteraceae 


3698 


102 


5 


8 


25 


3 


33 


14 


0.1% 


Proteobacteria 


Unresolved (Gamma-) 


10193 


8 


11 


1 


26 


79 


38 


17 


0.1% 


Proteobacteria 


Desulfovibrionaceae 


3473 


36 


12 


28 


6 


2 


8 


57 


0.1% 


Proteobacteria 


Endozoicimonaceae 


6801 


2 


7 


7 


1 


101 


5 


17 


0.0% 


Proteobacteria 


Rhodobacteraceae 


2007 


42 


6 


24 


9 


2 


13 


40 


0.0% 


Proteobacteria 


Endozoicimonaceae 


5020 


3 


1 


32 


6 


33 


49 


1 


0.0% 


Verrucomicrobia 


Verrucom icrobiaceae 


1910 


3 


3 


8 


6 


59 


20 


11 


0.0% 


Proteobacteria 


Rhodobacteraceae 


6387 


12 


8 


32 


14 


5 


20 


15 


0.0% 


Proteobacteria 


Endozoicimonaceae 


6815 


2 


1 


1 


30 


8 


43 


3 


0.0% 


Proteobacteria 


Unresolved (Gamma-) 


8480 


4 


3 


4 


2 


5 


3 


9 


0.0% 


Proteobacteria 


Unresolved (Alteromonadales) 



•Includes starved and unstarved animal samples. Numbers indicate the number of hits or abundance of each OTU (40,528 reads/sample). The first 9 OTUs in bold 
represent core OTUs that are part of the 1 2 most abundant OTUs described in Figure 1 . Full taxonomic lineage is defined for each OTU, where possible, in Supplemental 
Table 1. 

doi:1 0.1 371 /joumal.pone.0093386.t001 



the Ciona gut, as shared OTUs between distant populations are 
sufficient to mask the effects of their highly variable environments. 

The Ciona gut microbiota supports several familiar OTUs 

The observed differences in core OTUs between starved and 
unstarved Ciona samples (Figure 1; Figure SI; Table SI and S2) 
suggest that some core OTUs represent planktonic communities 
(e.g., Prochlorococcus) while others potentially form adherent 
communities in stable biofilms within the mucus layers (e.g., 
anaerobic Psychrilyobacter sp.). It is possible that adherent commu- 
nities benefit from early colonization of the gut surface at the onset 



of feeding and are influenced by host responses. Thus, coloniza- 
tion is likely influenced by physical properties of the Ciona gut (e.g., 
the presence of unique side chains in mucin-like glycoproteins 
[28]) as well as the types of bacteria colonizing planktonic or 
dietary particles. 

It is likely that some members of this microbiota, such as 
anaerobic Fusobacteria, stably associate with distinct compart- 
ments of the Ciona gut, and together with other members of the gut 
microbiota influence host fitness. Some surfaces of marine 
invertebrates can be colonized by specific groups of bacteria, such 
as members of the Oceanospirillales order [29-33]. In the Ciona 
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Table 2. Ciona core OTUs represent four phyla, 10 orders, and 13 families of bacteria. 





Core Phyla 


Order 


Family 


OTUs 


Bacteroidetes 


Flavobacteriales 


Flavobacteriaceae 


2 








2 


Fusobacteria 


Fusobacteriales 


Fusobacteriaceae 


2 


Proteobacteria 


Alteromonadales 


Pseudoalteromonadaceae 


1 






Psychromonadaceae 


1 


Shewanella 1 




Campylobacterales 


Campylobacteraceae 


2 




Desulfobacterales 


Desulfobulbaceae 


1 




Desulfovibronales 


Desulfovibrionaceae 


1 




Oceanospirillales 


Hahellaceae 


4 






Oceanospirillaceae 


1 




Rhodobacterales 


Rhodobacteraceae 


6 




Vibrionales 


Vibrionaceae 


4 








6 


Verrucomicrobia 


Verrucomicrobiales 


Rubritalea 


1 


35 



*OTUs that cannot be taxonomically sorted below class. 
doi:1 0.1 371 /journal.pone.0093386.t002 



gut, OTUs closely related (>97% identity) to Endozoicomonas sp., 
members of Oceanospirillales, can peak at 53.5% relative 
abundance (Table SI). Endozoicomonas sp. may serve essential roles 
in some marine invertebrates as they can dominate the surface 
mucus of many corals [31-35], be associated with internal tissue 
spaces [29], and be endosymbionts [36,37]. Diverse invertebrates 
appear to associate closely with these bacteria [30,38], which are 
predicted to serve important roles in metabolizing complex 
organic compounds [39] and even in sulfur recycling [40] , similar 
to other members of the Order Oceanospirillales. However, it 
remains unclear if they are found as part of other marine 
invertebrate microbiomes [3,41,42]. 

Overall, the abundance of core OTUs in Ciona can vary 
considerably. The abundance of a Psychrilyobacter (Fusobacteria)- 
related OTU, for example, can reach 30% in some samples. In 
addition, the abundance of several other OTUs within Bacter- 
oidetes and Gammaproteobacteria can vary in the Ciona gut even 
among starved animals (Table SI). Bacterial communities that 
associate with some coastal invertebrates are subject to anthropo- 
morphic influences [31,32]. We surveyed marine invertebrate- 
associated microbial communities by screening 16S sequences in 
the available databases and recorded OTU matches (>97% 
identity) in various animal taxa (Table S3). More specifically, Ciona 
core OTUs were compared to coral gastric cavity [43], polychaete 
intestine [42], abalone intestine, and two sea cucumber intestinal 
microbiomes [44,45] and the sharing of some bacterial families 
and genera were noted (Table 3). However, this sharing was not 
across any significant number of genera, with the exception of one 
sea cucumber, which shared 12 OTUs among 5 unique genera 
with Ciona (Table 3). About 29% of the Ciona core OTUs (10 of 35) 
have been identified in the gut of various sea cucumbers species 
(although not all in any one species; Table S3). Indeed, recent 
work has revealed that many invertebrates maintain bacterial 
communities of related (upper-level) taxonomic groups composed 
of shared classes and families within Proteobacteria, Fusobacteria, 



Bacteroidetes, Actinobacteria, and Firmicutes in their guts. 
Whereas a few of the Ciona core OTUs previously have only 
been described in environmental samples (Table S3), most of the 
core OTUs have been described as associated with other animal 
taxa, for example, the Fusobacteria genera, Psychrilyobacter and 
Propionigenium. Taken together, these findings suggest that some 
bacterial families form abundant communities in the gut 
environment of diverse marine invertebrates, possibly demonstrat- 
ing a preference for the feeding habits and/ or the gut ecosystem of 
these animals. 

Illumina sequencing reaffirms OTU findings from clone 
library surveys 

Although the power of the studies reported here is somewhat 
limited by the sample size, it is noted that between 2010-12, 16S 
rRNA PCR clone library surveys of 70 animals (46 from WH, 20 
from SD, and 4 from Naples) also were conducted. Samples were 
collected randomly throughout each year; some were starved, 
others had full guts and others had guts devoid of fecal matter. 
Approximately 3400 Sanger-based sequencing reads (nearly 1 700 
clones sequenced in both directions) resulted in 1617 near full- 
length 16S sequences (Genbank accession KF798359-KF799975) 
resulting in 442 unique OTUs. The 47 most abundant OTUs 
(seen more than >5 times) are markedly similar to the core 
community that is predicted by Illumina sequencing of the starved 
and unstarved sample sets (Figure 2; Table S4; Table S5). This 
finding reaffirms that the core bacterial community within the 
Ciona gut is stably associated, abundant, and can be captured 
irrespective of handling or sampling method. 

In interpreting host-OTU relationships between collection sites, 
it should be noted that recent evidence indicates Ciona intestinalis 
could be experiencing a speciation event currently defined by two 
subtypes [46-49] . The ' A' subtype is represented by populations in 
the Mediterranean and the west coast of North America (here 
Naples and SD, respectively); the 'B' subtype is found along 
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Table 3. Comparison of the 35 Ciona core gut OTUs to other marine invertebrate gut microbiota. 





Subject organism 1 


Gut OTUs 
reported 


Phyla 
reported 2 


Phylun 
hits 3 


i-level Family-level 
hits 3 


Genus-level 
hits 3 


Unique Genus 
hits 


Coral {Galaxea fascicularis) gastric cavity 


17 


3 


11 


6 


4 


1 a 


Polychaete (Neanthes glandicincta) intestine 


84 


6 


32 


15 


5 


1 b 


Abalone {Haliotis discuss hannai) intestine 4 


63 


5 


56 


50 


36 


2 C 


Sea cucumber (Holothuria leucospilota) intestine 


139 


3 


53 


48 


48 


2 d 


Sea cucumber (Apostichopus japonicus) intestine 


217 


8 


131 


110 


12 


5 e 



All comparisons are between Ciona gut core OTUs and other marine invertebrates. Note that priming strategies and sequencing methods vary greatly. Therefore, OTU 
similarities are approximate however at >95% confidence over the regions compared, dumber of OTUs sharing taxonomic hits (with Ciona core OTUs) in BLAST 
database searches. References cited in the text. Unclassified bacteria not shown. 3 Number of subject OTUs shared with Ciona core OTUs at each taxonomic level 
specified. 4 Unpublished data deposited in NCBI. a Endozoicomonas; b Vibrio; c Vibrio and Shewanella; d Shewanella and Pseudoalteromonas; e Vibrio, Rubritalea, Arcobacter, 
Pseudoalteromonas, and Shewanella. 
doi:1 0.1 371 /journal.pone.0093386.t003 



Scandinavian coastlines and along the northeastern coast of North 
America (represented here by WH). Both subtypes are geograph- 
ically restricted and demonstrate distinct ranges of temperature 
tolerance. Recently, a unique assortment of microbiota inhabiting 



the tunic of Ciona intestinalis was described from 'B' subtype animals 
(WH samples) [50]. Notably, almost no overlap (i.e., shared 
OTUs) exists with the microbiota of the gut. 




Figure 2. Unrooted distance tree comparing lllumina-derived core OTUs and abundant Sanger-sequenced clone library 16S 
samples. Pairwise comparisons of OTUs indicates that a significant proportion of the core OTUs were captured consistently over a three year 
collection period spanning all seasons. 
doi:1 0.1 371 /journal.pone.0093386.g002 
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Selection of a core microbiota may not be limited to 
feeding habits and environment 

A complex assemblage of a stably associated core bacterial 
community has been defined in the gut of Ciona intestinalis, sampled 
from three disparate geographic locations. To our knowledge, this 
is the first report of a core gut microbiota from a marine 
invertebrate under starved and unstarved conditions. Ciona may 
have evolved some level of symbiotic interdependences, possibly 
extending the range of host physiological responses, which could 
be targeted to identify symbiotic interactions across the gut 
epithelium. Studies of the host-microbial dialogue of the gut could 
help reveal basic requirements that drive and/or help sustain 
immune homeostasis of the gut. Immune-type molecules, such as 
the immunoglobulin variable region-containing chitin binding 
proteins (VCBPs), are expressed in high abundance in the gut 
epithelium and it remains to be determined what role, if any, these 
molecules play in modulating the composition of closely-associated 
gut microbial communities. VCBPs are secreted proteins that have 
been shown to bind bacteria both in vivo and in vitro; because they 
associate closely with epithelial processes and microvilli in the gut 
lumen [15], it is most likely that they interact with bacteria 
adhering to the gut mucosa. Determining the precise mechanisms 
by which Ciona establishes and maintains its gut microbial 
ecosystem, from both anatomical and immunological perspectives, 
could have important implications for our understanding of the 
evolution of gut immunity and the nature of innate immunity at 
the gut epithelial interface [14]. 

Supporting Information 

Figure SI Summary of the number of core OTUs shared 
between all (left), not starved (center), and starved 
(right) samples. Bins are the number of samples necessary for 
an OTU to be detected and be considered "core" and the height is 
the number of OTUs in that bin. 
(PDF) 
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